Vortex characterisation of frustration in the 2d Ising spin glass 
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The frustrated Ising model on a two-dimensional lattice with open boundary conditions is revis- 
ited. A hidden IZ2 gauge symmetry relates models with different frustrations which, however, share 
the same partition function. By means of a duality transformation, it is shown that the partition 
function only depends on the distribution of gauge invariant vortices on the lattice. We finally show 
that the exact ground state energy can be calculated in polynomial time using Edmonds' algorithm. 
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Introduction Spin glasses [l[ are magnetic materials 
in which the magnetic moments are subject to ferro- 
magnetic or anti-ferromagnetic interactions, depending 
on the position of the moments within the sample. The 
system is frustrated in the sense that the arrangement of 
spins which minimises the total energy cannot be deter- 
mined by considering a local set of spins. Stated differ- 
ently, the change of a single spin might cause a reorder- 
ing of many spins when the system relaxes towards a new 
minimum of energy Q. Spin glasses undergo a freezing 
transition to a state where the order is represented by 
clusters of spins with mixed orientations. The relaxation 
times towards equilibrium are typically very large, which 
impedes efficient simulations. 

Many efforts have been undertaken to explore equilib- 
rium properties of spin glasses by means of Monte Carlo 
simulations 

(MS- Thereby, many insights have been 
obtained from the simple case of the 2d Ising model on 
a square lattice. For the discrete model, the bond in- 
teractions take values ±1 at random, and the model is 
characterised by the the probability re of finding an anti- 
ferromagnetic interaction at a given bond. 

As first noticed by Bieche et al. [ll| and further elabo- 
rated by Nishimori |l2l . [l3j , the Ising model with random 
distribution of anti-ferromagnetic bonds has a hidden 
gauge symmetry. As discussed below, this symmetry im- 
plies that gauge invariant observables such as the thermal 
energy or the specific heat are unchanged by a certain re- 
distribution of the anti-ferromagnetic bonds (which may 
also change considerably in number). By exploiting this 
invariance, Nishimori was able to obtain exact results for 
special distributions of the bond frustration and temper- 
ature (HB3. 

In this letter, we further explore the consequences of 
the hidden TLi gauge symmetry which relates models with 
the same partition function but with different frustra- 
tions. We show that within the class of gauge equivalent 
models, there is always one model for which the ground 



state is homogeneous. Furthermore, we show by means 
of a duality transformation that the partition function of 
any spin glass only depends on the distribution of gauge 
invariant vortices. The amount of frustration in the 2d 
Ising model is quantified by counting the fraction p of 
vortices (non-trivial plaqucttes) in a given bond distri- 
bution, and we determine the exact ground state energy 
as a function of p using Edmonds' algorithm [Til Il5|. 

Hidden gauge symmetry The partition function of the 
frustrated Ising model involves a summation over all spin 
configurations {cr x } 



Z = 



(1) 



e={xy) 



where the spins located at the sites of the lattice take 
the values o~ x = ±1. We work with open boundary con- 
ditions throughout this paper. The sum in the exponent 
extends over all bonds i = (xy) and the bond variables 
Ul are chosen equal to +1, except for a fraction re of 
randomly distributed bonds where Ut = (—1). In the 
zero temperature limit, /3 — > 00, the anti-ferromagnetic 
couplings induce frustration. 

It was first observed by Nishimori [12, HH that bond 
distributions with vastly different values for re may still 
share the same thermodynamical properties. This is due 
to a 7Li gauge symmetry. The partition function, eq. ([IJ , 
and observables such as the thermal energy, 



E(/3) 



d\nZ 

d[3 



i=(xy) 



(2) 



are invariant under the following change of bonds and 
spin variables: 



U (xy) 



Q,(x) o~(x) 

n(x) U( xy ) n(y) 



(3) 



where the gauge transformation takes values in W12, 
Q(x) = ±1. 
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A simple consequence of this symmetry is that we can 
always find a gauge for which the ground state is uniform. 
In particle physics, this gauge is known as Landau gauge. 
It corresponds to the bond distribution with the minimal 
number of anti-ferromagnetic bonds: 



E 

i=(xy) 



a n 



(4) 



For this choice of bond distribution, the uniform state is 
a ground state: 

— <jy L — const. \fx,y, (Landau gauge) . (5) 

To see this, we can use eqs. (J3J and ((3]) to express the 
energy of a given spin configuration {a} in a Landau 
gauge background as 

E[a] = Y, a * U ? L a v = ~ E U ^ L 

l=(xy) (=(xy) 

> - E U ? L = E & U ? L o^=E[o a -] 



l={xy) 



£={xy) 



where we used the maximum condition eq. ^ for the in- 
equality and the definition cq. (JS|) for the uniform ground 
state, which satisfies a^ L a^ L = 1, Vx, y. It might turn 
out that several sets of gauge transformations {H,} ac- 
complish the (global) maximum in The degeneracy 
of the spin glass ground state is then twice the number 
of these sets. What we have shown is that among these 
sets there is always the uniform state. 

Vortex description Let us further specify the frustra- 
tion of the model in a gauge invariant way. To this end, 
we introduce plaquette variables, U p = Y\i^ P Ui, con- 
structed from the given bond background. (The product 
is along the four bonds forming an elementary square 
p.) This definition is borrowed from lattice gauge the- 
ory, where a non-trivial value U p = — 1 indicates that a 
vortex intersects the elementary square p. Formally, 
we set Vp* = Up for all dual plaquettes p* of dimension 
(d — 2) and locate a vortex if V p * = — 1. (In d = 2, 
p* is the midpoint of a square p.) With open boundary 
conditions the plaquettes or vortices determine the bond 
variables Ut completely, up to TLi gauge transformations. 
For this reason, the distribution of vortices is the proper 
measure to unambiguously quantify the frustration of the 
model. 

To see this in more detail, we invoke a duality trans- 
formation. We begin by expanding the exponential in (fTJ) 
and use [a x Uia y ] 2 = 1 to find 



(cosh/3) 



N, 



E 

{<j x } t={xy) 



[l + tanh(/3) a x U t a v ] , (6) 



where Ni is the number of bonds. Expanding the prod- 
ucts and summing over all configurations of TV spins as 
usual, the partition function becomes a sum over all sets 
of closed loops on the lattice: 



Z 



2 w (cosh/3) Arf Y [tanh/3] L ^ W(tf) 



(7) 



Notice that a particular loop set in the sum may, in 
general, contain multiple closed loops that are discon- 
nected, touching or even intersecting, i.e. the loops in 
the set may share points, but not bonds. Furthermore, 
L(^) = |^| is the total number of bonds in the loop set, 
while W(^), which is called the Wilson loop in gauge 
field theories, is their product: 



W{ff) = ]Ju t 



(8) 



The loop formulation eq. (J7) is valid in any number 
d of space dimensions. It is remarkable that the bond 
variables Ut of the given spin glass instance only appear 
in the gauge invariant Wilson loop factor W^). To show 
that this information is equivalent to the distribution of 
non-trivial plaquettes, we only have to appeal to Stokes' 
theorem 



W(V) 



p j 



u p = ]Ju e . 



(9) 



where A is any area on the lattice bounded by the loops 
in c 6 '. Thus, we have shown that the distribution of non- 
trivial plaquettes U p (or vortices V p * ) determines the Wil- 
son factor for all loop sets and thus characterises the 
given spin glass completely via eq. 0. 

To complete the gauge invariant duality transforma- 
tion, we would like to convert the loop formulation ([7]) 
to a spin model on the dual lattice, with the vortices 
appearing as the spin glass background. More precisely, 
we map loop sets onto configurations of spin variables 
t p * <S {±1} as follows: 

1. start with an empty configuration r p * =1 Vp* ; 

2. decompose ^ in a union of connected loops ^ ; 

3. for each choose an area Ai^i) bounded by ^ ; 

4. transform the spin configuration {t p * } according to 



Tp* 



-t p , ; p e A(%) 
p £ A{%) ; 



(10) 



5. repeat steps 3. and 4. until all ^ have been used. 

The resulting spin configuration {r p * } does not depend 
on the choice of loop decomposition in step 2, but it 
does depend on the subsequent choice of surfaces A(%) 
in step 3. An exception is the case d = 2 where Ai^€i) is 
the unique area enclosed by the loop ^ and the resulting 
{r p » } is unique. On the other hand, the reconstruction of 
the ^ from the dual spin configuration is unique in any 
number of dimensions d, since {r p * } has the property 



n 

p*£dt 



V 



-1 ; ietf 
-1 : H c €. 



(11) 
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This is because n p * 



\ [l — r p »] indicates whether a 
particular plaquette p is contained an even (n p * = 0) 
or odd [n p - = 1) number of times within the total area 
A(^). A bond I on the initial lattice is clearly part of 
the boundary ^ iff the "area count" of all plaquettes 
attached to it is odd, 

p*£dl* p*£dt* 

from which cq. (fTT]) follows using (— l)^ 1-2 ) = z, z g Sa- 
Thus, the map — > {r p *} is one-to-many except for 
d = 2, where it is one-to-one. 

With the construction of the dual spin variables, the 
duality transformation is now simple. First we have (in 
any number of dimensions), 

*W = Evr>=E Jt^w] (12) 



(p*9*> 



<P*9*> 



= u p=i[ u p f * = n u P li ' Tpt] ■ 

peACtf) p p 

In d = 2, the map ^ — > {Tp*} is one-to-one, and we can 
directly replace the sum over loop sets by a sum over dual 
spin configurations: 



2 Ar -^ 2 sinh(2/3) A ^ 2 x 



(13) 



exp 



<P*9*> 



where /3 = — In (tanh /?) /2 and Vp* = U p is the vortex 
background. This formulation lives entirely on the dual 
plaquettes, i.e. the sites of the dual lattice. For conve- 
nience, we can also rewrite the Wilson loop factor as 



n 



(14) 



p : 



:V„.=-I 



This shows explicitly that the partition function of the 
spin glass is entirely determined by the location of vor- 



tices (V p - 



T) on the dual lattice. 



Finally, we can also rewrite the vortex contribution to 
the partition function as an imaginary external magnetic 
field. However, although the low-temperature limit is 
attained for (3 — > suggesting a standard strong-coupling 
expansion, the dual theory is still tedious to solve due to 
the inherent sign problem. 

The extension to d > 2 space dimensions is straight- 
forward, with spin configurations now assigned to dual 
bonds, plaquettes etc. However, the transformation (JTUJ) 
is complicated by the fact that a given loop configura- 
tion ^ does not lead to a unique dual spin configura- 
tion { t p * }, because there are many surfaces A(^) over 
^ which enter the construction of the dual spins. In a 



certain sense, this is a new (geometrical) gauge symme- 
try which requires a unique prescription of how to assign 
areae A( c £) to closed loops in d > 2. Only such a fixed 
partition function is equivalent to the loop formulation 
cq. ©, but the "area fixing" is presumably NP hard in 
d > 2. 

Numerical studies In d = 2, the ground state energy 
can be obtained in polynomial time. To see this, consider 
the dual formulation (|1 3|) in the low-temperature limit 
P — > 0, when the exponential can be expanded. Using 
the representation ([14]) for the Wilson factor, the subse- 
quent sum over all spin configuration gives zero unless 
every spin variable r p » appears an even number of times, 
i.e. all spin variables from eq. (|14p have been paired by 
corresponding variables from the expansion of the ex- 
ponential. At order , the exponential gives one pair 
of adjacent r p * 's which can match two adjacent vortices 
from the Wilson loop; at order f3 2 , we can match two 
pairs of adjacent vortices, or one vortex pair separated 
by two dual bonds. Proceeding in this way, it becomes 
clear that the lowest non-vanishing contribution to the 
partition function is given by the minimal matching of 
all vortices in the Wilson factor, 



J d=2 



,JV- 



? / 2 sinh(2/3)" f / 2 -v-P Na 



where Na is the length (number of bonds) in the mini- 
mal matching, and v is the degeneracy (number of dis- 
tinct minimal matchings). The ground state energy thus 
becomes 

E = -N e - N A lim = 2N A ~ N e . (15) 

The core of the solution is therefore the minimal match- 
ing of vortex pairs on the dual lattice, for which Ed- 
monds ' algorithm [Til Il5j gives an answer in polynomial 
time [ll|]. 

Figure Q] shows a random distribution of Ny = 150 
vortices on a L x L lattice with L — 30 and open bound- 
ary conditions, and the corresponding minimal number 
of frustrated bonds obtained with Edmonds' algorithm. 
Also shown is the ground state energy per bond as a 
function of the vortex density p = Ny/(L — l) 2 . The 
data comprise an average over 100 random vortex distri- 
butions for each value of p. It turns out that the ground 
state energy per bond is well fitted for small values of p 
by 



Eo( P )/N e 



-1 



0.39(2) p 0A9 ^ 



(16) 



Two limits can be obtained analytically: For p = 0, we 
recover the ferromagnet with Eo(0)/Ne = — 1. For the 
other extreme p = +1 and for L odd, each plaquette of 
the lattice carries a vortex. Each vortex pair is saturated 
by one frustrated bond. Hence, Na = (L — l) 2 /2 bonds 
out of 2L(L — 1) bonds are frustrated so that 



Eo(l)/N e 



2N A - N e 



1 

2L 



;i7) 



FIG. 1: Left panel: Ny = 150 randomly distributed vortices (red squares) on a 30 x 30 lattice and the corresponding minimal 
matching with anti-ferromagnetic bonds (black lines). Eq/Ni = —0.82759(1). Right panel: Ground state energy per bond as 
a function of the density of frustration p on 60 2 , 90 2 and 120 2 lattices with open boundary conditions. The two fits refer to 
eqs. (|16l) and (jTHJ) , respectively. 



Here we obtain a 1/L correction to the infinite volume 
result. Finally, (fLo)) implies that the derivative dE /dp is 
singular for vanishing vortex density indicating that per- 
turbation theory with respect to the frustration density 
is not always justified. 

There is a simple way to understand cq. (|16[) . Suppose 
that the Ny vortices are distributed uniformly (without 
correlations) on the lattice. Each vortex occupies, on 
average, an area which should be taken (because of 
the taxi driver metric) as a rectangle with side length 
p -1 / 2 . The greedy matching of two neighbouring vortices 
then has length p^ 1 ^ 2 , and the complete greedy matching 
of all vortices requires Na — Nv /2 ■ p^ 1 ^ 2 bonds. Hence 
we obtain from eq. (|15[) a ground state energy of 



"aT" 1 + VP 2L(L - 1) 



1 



+ (18) 



This explains the peculiar ^fp— behaviour for uniformly 
distributed vortices. The true ground state is, of course, 
based on the minimal, not the greedy matching, i.e. the 
prefactor ^ji in the above estimate must be multiplied by 
the ratio a < 1 of the bond numbers in the minimal and 
greedy matching. At small vortex densities p C 1, the 
empirical law (fTB")) suggest a « 0.8, while the exact result 
(fl7| implies a = 1 at p = 1. 

Both formulae eq. (fTB")) and eq. (jT5J) give good estimates 
of the ground state energy for spin glasses with a uniform 



(uncorrelated) vortex distribution. To test an extremely 
correlated situation, we place Ny vortices as Ny/2 pairs 
of exactly one bond distance (dumbells). Obviously, the 
greedy matching is also minimal in this case, i.e. we have 
one matching bond per vortex pair, or Na = Ny/2. The 
exact ground state energy per link is thus Eq /Ni = — 1 + 
Ny/Ni = — 1 + p/2 for L 3> 1. By comparision, the 
estimate (|18f has a maximal deviation of 14.6% for this 
configuration, attained at p ps 0.34, while the fit eq. (fTB")) 
intended for small p has a relative error of less than 9% for 
all p < 0.6 even on the correlated dumbell configuration. 
Thus, it seems that the above estimates work reasonably 
well even for correlated vortex distributions. 

Conclusions In conclusion, we have emphasised the 
importance of a gauge invariant classification of frustra- 
tion. The partition function was found to depend solely 
on the distribution of gauge invariant vortices on the lat- 
tice. In two dimensions, the exact ground state energy 
was calculated in polynomial time using Edmonds' al- 
gorithm. Using numerical simulations, we find the de- 
viation of the ground state energy from that in the fer- 
romagnetic case is to a good extent proportional to the 
square root of the vortex density. 
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